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Abstract 

We study the oscillation frequencies and radiative decay rates of surface plasmon modes of a simple-metal sphere 
as a function of sphere radius without any assumptions concerning the sphere size. We re-examine within the 
framework of classical electrodynamics the usual expectations for multipolar plasmon frequency in the so called 
"low radius limit" of the classical picture. 
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1. Introduction 

The dielectric properties of metals, as well as 
those of semiconductors with high electron concen- 
tration, are due to collective effects arising from 
the Coulomb interaction between charges. In sim- 
ple bulk metals the conduction electrons can be 
considered as a free-electron plasma. Frequency de- 
pendence of some of optical properties can be well 
described at a quantitative level by the Drude- 
Lorentz dielectric function [1] . Optical properties 
of the electron gas in bulk metals, in proximity 
semi-infinitive surfaces, in thin films, and in metal- 
lic particles can be characterized by the eigenfre- 
quencies of the system depending on free electron 
density and the geometry of the system. If we talk 
of " plasmons" or " plasma waves" , we mean eigen- 
modes of the self-consistent Maxwell equations for 
the system in the absence of an external electro- 
magnetic field (or in a direction orthogonal to the 



field) (e.g. [2], [3]). " Surface plasmons" , are used as 
a name for electromagnetic eigenmodes which are 
maximal near the surface. The time dependence 
of eigenmodes of a free-electron system is charac- 
terized by corresponding eigenfrequencies with the 
real part defining the frequency of oscillation, and 
the imaginary part defining the radiative damping. 
Usually the eigenmode problem of a metallic 
sphere is studied in the limit of very small size pa- 
rameter (retardation effect omitted) , e.g. [4], [3], [5] 
and the radiative damping of plasmon oscillation 
is not included. The dipole mode eigenfrequency 
is then expected to be equal to cjp/\/3 being re- 
sponsible for the "giant dipole resonance" result- 
ing from the Mie scattering theory ([6], and also 



e.g. [7], [4], [8], 



is the modified plasma fre- 



quency, which can include or not the cluster core 
polarizability and (or) the spill-out effect of elec- 
tron density at particle border, depending on the 
model approximations in effect. 
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In [10], [11] we have reconsidered the eigen- 
value problem of a free-electron metal sphere as 
a function of sphere radius without any assump- 
tion concerning the particle size and including 
higher eigenmodes than the dipole ones. We have 
studied the dipole (1 = 1) and the higher polar- 
ity plasmon eigenfrequencies uji{R) as well as the 
plasmon radiative decay rates uj'/lR) as a func- 
tion of the particle radius R with no assumption 
concerning the lower limit of the particle size in 
numerical modelling (retardation effects included) 
for I = 1,2, ...6. In [11] we have also studied the 
plasmon manifestation in scattering and absorb- 
ing properties of the sphere of arbitrarily large 
size (retardation included) within full scattering 
Mie theory. 

In the present paper, we use the same "exact" 
solutions of the eigenmode problem for I ^ 1,2, ...6 
and 7, 8... 10 in addition, and re-examine the usual 
expectation for multipolar plasmon frequencies in 
the so called " low radius limit" . If the particle is 
formed from ideal metal (free electrons do not suf- 
fer from collisions 7 = 0) and is embedded in vac- 
uum {eout = 1) the multipolar plasmon frequen- 
cies according to the "low radius limit" approxi- 
mation are expected to be ujqj ~ uipy^l/ {21 + 1) 
(e.g. [12], [13], [3], [5] ). The well known dipole mode 
frequency ujp/\/3 is obtained for 1 = 1, while for 
increasing / the eigenmode frequencies approach 
the frequency of plane surface plasmon at u!p/\/2, 
in spite of the fact they result from the "low ra- 
dius approximation" (i.e. from the limit of i? ^ 0, 
while plane surface limit is i? ^ oo). In this paper 
we study the reasons of underlying causes for this 
paradox. 



2. Formulation of the eigenvalue problem 
for a sphere of arbitrary size 

The starting point is provided by the self- 
consistent Maxwell equations: 



VxE=-f 



V • E = p/eo 
VB = 



(1) 



with no external sources: pe^t — 0, jext = so 
j and p are induced current and charge densities 
respectively. 

The frequency dependent dielectric function 
s(lo) = Si„{u)) and conductivity a{uj) = am (to) of 
the sphere is assumed to have the constant bulk 
value up to the sphere border. The dynamic, lin- 
ear response of the sphere material is described 
within standard optics, so the local proportion- 
ality between the electric displacement D and 
electric field intensity E at the same point in 
space are valid: D(r, oj) =E(r,a;) -I- — j(r,u;) = 

eo(l + ^f^)E(r,w)=eoe(a))E(r,cj). The sphere 
is embedded in nonconducting and nonmagnetic 
medium a out = and s{uj) = Sout will be as- 
sumed to be Eout = 1 in all numerical illustra- 
tions. The dielectric function of the sphere will 
be assumed to be the Drude dielectric function 



.H = i- 



-t^UJ 



-. We look for solutions fulfill- 



ing Maxwell's equations in the form of transversal 
waves (V • E = 0) in two homogeneous regions 
inside and outside the sphere so the wave equa- 
tion: V2E(r) + V(V • E(r)) - ^^ = for 
harmonic fields E(r,cj) = ^e^^C^r-^^t) j-gj^^gg ^q 
the Helmholtz equation: 



V2E(r)+g2E(r) = 



(2) 



where: q = qin inside the sphere, q = qout in 
the sphere surroundings, q^ = qoy/si^ , qout = 
qoy/Sout and qo = —■ The well known scalar solu- 
tion of the corresponding scalar equation (e.g. [4], 
[8]) in spherical coordinates (r, 9, (/>) reads: 



ipim{r, 0, (/)) = Zi{qr)Yirn{0, (p). 



(3) 



where I = 1,2,..., m = G,±l,...,±/, Yim{9,(t>) 
are spherical harmonics, and Zi{qr) are spherical 
Bessel functions jiiqinr) inside the sphere and the 
spherical Hankel functions hi{qoutT) outside the 
sphere. 

Because various notations have been employed 
in different papers and textbooks and none appears 
to have general acceptance, let's recall that the 
spherical Bessel functions: ii{z) = y^^J;+i(z), 

and hi {z) = ji {z)-i- ni{z) = Jj^H^Ji (z) where 

^li-^) ~ \f^-^i+-i^)- ^^'^ functions J;^i(z), 



H^ i{z) and Ni^i{z) are Bessel, Hankel and 



Neuman cylindrical functions of half order of the 
standard type according to the convention used 
e.g. in [7]. 

From scalar solution ipim one can construct two 
independent solutions of the vectorial wave equa- 
tion (2), one with vanishing radial component of 
the magnetic field: 



E(r) = Birn{l/q)V X V X [riPim), 
H(r) = Bi,n{q/iqo)^ x (rV'/m). 



(4) 
(5) 



and the other with vanishing radial component 
of the electric field: 



E(r)=A,„V X (r^pim), 

H(r) = Ai^{l/qo)V X V X (r^im), 



(6) 
(7) 



Aim and Bim are constants that take different 
values A^;^ and B;™ inside and Af^* and BZ* out- 
side the sphere. The explicit expressions for the 
solution with the nonzero radial component of the 
electric field Er y^ (and the magnetic field tan- 
gent to the sphere surface Hr — 0), which is named 
transverse magnetic (TM) mode in analogy to the 
fiat surface interface case {p polarization, or " elec- 
tric wave" in terminology of [7]) read: 



Er{r,0 
Hg{r,9 



^Biml{l + l){qr)-^Zi{qr)Yim{eA), 
^BUqr)-'[qrZi{qr)]'dYim/d0, 
^Bimim(qrsiTLey^[qrZi{qr)]'Yim(0,(t)), 
= 0, (8) 

^ Bim[e{uj)Y'^{m/ sme)Zi{qr)Yim{e,^), 
^iBim[e{uj)Y'^Zi{qr)dYimlde, 



The expression for the orthogonal solution with 
-Er = results from eqs.(6,7) (and is named trans- 
verse electric (TE) mode in analogy to the fiat sur- 
face interface case {s polarization)). The prime in- 
dicates differentiation in respect to the argument, 
which is qinr or qout^ correspondingly. We focus 
our attention on TM mode only. 

The continuity relations at the sphere boundary 
for the tangential components of the electric field 
(the continuity of Eg and E^p) lead to the same 
condition: 



while the tangential components of the magnetic 
field (the continuity of Hg and H^) lead to the 
condition: 

BZ^^3i{zb) - BZ'V^:^Mzh) (10) 

where: 



ZB — qinR — —Ry/^i 



(11) 



is the argument of the Besscl function ji , and 



zh = qoutR 



— ti^Eout — Zb — J— 



(12) 



is the argument of the Hankcl function for r ~ 
R. The continuity relations for TM mode lead to 
non-trivial solutions (e.g. non-zero field amplitudes 
Elm inside and outside the sphere) only when: 



[zBJiizB)]' _ [zHhiizn)]' 

£injl{zB) £outhl{zH) 



(13) 



We are interested in the properties of the sphere 
in the frequency regime of anomalous dispersion 
£m(w) < 0. In that region only the TM eigenmodes 
exist, while the equation dispersion relation for TE 
mode has no solution for ei„(w) < ([5] or [14]). 
Zi{qr) = jiiqinf) is then a function of a complex 
argument and the solutions given by eqs. (8) are 
called "surface modes". The fields are maximal at 
the sphere surface, with exception of the / — 1 
mode which is uniform throughout the sphere ([4] 
or [14]). 

On writing down the dispersion relation for 
the TM mode (13) in terms of the more com- 
pact Riccati-Bessel function ipi (z) ~ z ■ ji{z) and 
^i (z) = z ■ h\ (z), the dispersion relation for the 
TM mode reads: 



/£outS.l [zh) ^'l (zb) 



r^,(zB)c;(2ff) = o(i4) 



BZ{zb)-'[zbJi{zb)]' 



BZ\zH)-^[zHhi{zH)]'{9) 



The boundary conditions are then satisfied 
only by a discrete set of characteristic complex 
values zi which are the roots of the complex 
fu nction Di{z) = ^SouS {zHi(^))i^'i (zBi^j)) - 
sj£in{u:)i)i {zB{i^))i!i {zh{uj)) of complcx argu- 
ment z — z{ijj,R). Discretization of complex roots 
zi means the discretization of corresponding val- 
ues uj = $7;, I = 1,2,3... which are allowed to be 
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Fig. 1. Plasrnon oscillation frequencies iOi{R) as a function 
of sodium sphere radius R for Z = 1, 2, ...10 (rigorous solu- 
tion). The picture illustrates the coincidence of the plas- 
rnon frequencies LOi{Rmin,l) with the corresponding value 
a;o I obtained within vanishing size approximation (open 
circles). 7 = 0. 

complex: fli = LUi+iuji . They define discrete eigen- 
niode frequencies uji and damping rates w" for the 
TM mode being the sum of corresponding compo- 
nents of (8) multipHed by e'"'* = e''^'*e""*. The 
analytic form of z; = zi{Cli{R), R) is not known, 
nor the analytic form of the relation ni{R). Let's 
notice, that neither znit^) nor zb{lo) separately 
are appropriate to define the set of discrete char- 
acteristic values, contrary to what is suggested in 
[8]. 

We solved the dispersion relation (14) with re- 
spect to il; numerically by treating the radius R 
as an external parameter. Riccati-Bessel functions 
^i, xi and ^i (and their derivatives with respect to 
the corresponding arguments zh and zb) were cal- 
culated exactly with use of the recurrence relation. 

We have used the Mueller method of secants of 
finding numerical solutions of the function f{v) = 
when one knows the starting approximated val- 
ues lying in the vicinity of the exact function pa- 
rameter w, which can be complex (the "root" func- 
tion of the Mathcad program). For given I and 
given R, the complex eigenvalue Qi was treated 
as the parameter to find, successive values of R 
were external parameters and where changed with 
the step Ai? « 2nm up to the final radius value 
R = 300nm. The values for uJi{R) and u>i{R) were 



Fig. 2. Rate of plasmon radiative damping lu',' (R) as a func- 
tion of sodium sphere radius R for / = 1,2, ...10 resulting 
from non-approximated radius dependence for 7 = 0. 

searched for by starting from approximate values 
of the root procedure chosen from the range from 
WpvS up to ujp\/2 correspondingly and for nega- 
tive values of Wj". The numerical illustrations have 
been made for a sodium sphere described by the 
Drude dielectric function with LUr, — 5.6 eV. 



3. Results 

Very careful study of roots of the function 
Di{Vti) of parameter fli{R) for given I for the de- 
creasing limit of radii R leads to the conclusion, 
that if the sphere is of the radius smaller than the 
characteristic radius Rmin,h there exist no ^i{R) 
real nor complex. So the complex eigenfrequen- 
cies fti{R) = u!i{R) + iuji{R) can be attributed to 
the sphere starting from the characteristic radius 
R = Rmin,i 7^ in given /. There exist no purely 
real solution for Qi: surface plasmons are always 
damped, even if the dielectric function s{uj) is real 
(7 = 0). 

Figure 1 and 2 (solid lines with closed spheres) 
illustrate the obtained loi{R) and uj'i'iR) depen- 
dencies for 7 = and / = 1, 2, 3, ...10 starting from 
Wi(i?min,i) and to'i'iRmin.i) valucs. These figures 
complete the picture for the R -^ Rmin,i limit of 
the corresponding dependence presented in [10], 
[11] for / = 1, 2, ...6, figures 1 and 3, while in [10] we 



4,0 



> 

o 3,6 



S"" 



3,4 



3,2 



3,0 



1 = 9 1=10 



1 = 2 



:i=i 



• Y=OeV 
■ Y= 1 eV 



20 



40 

R 



60 

[nml 



80 



100 



Fig. 3. The relation of the plasmon frequency uji and -Rmin,; 
for successive values oil = 1, 2, ...10 for electron relaxation 
rates 7 = and 7 = 1 eV 




Fig. 4. The dependence of damping rates: u>"{R) for the 
electron relaxation rate 7 = 1 eV for successive values of 
I = 1,2, ...10. 

did not study the limiting case of uji{R -^ Rmin,i) 
nor oj'i'iR -^ i?min,0 in detail. More careful search 
for these frequencies in the limit of smallest sphere 
still characterized by the eigenvalues w/(i?min,i) 
have shown, that they tend to the values which 
can be approximated by loqj values: 



^liRmin,l) ~ ^0,1 — ^f 



I 



21 + 1' 



(15) 



as illustrated by the hollow circles in figure 
1. Our numerical experiment shows that -Rmin,; 
dependence on / can be described as i?min,i ~ 
C [I {21 + 1)] ' with the proportionality constant 



C depending on density of free electrons. Rmm,i 
can be e.g.: i?min.(=4 = 6nm, but it can be as 
large as i?min.i=io = 87.2nm (the size parameter 
2ttR/X ~ 1 for optical wavelength A). 

The frequencies loq.i result from the dispersion 
relation (14) in the limit of small size parameter of 
the power series expansion of the spherical Bessel 
and Hankcl functions. 



hi{z) 



z' 


1 


(2Z + 1)!! 
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(16) 
(17) 



where {21 ± 1)!! = 1 x 3 x 5 x ... x (2? ± 1). If 
one employs the widely used rough approximation 

(e.g.[3], [5], [14]): 



'ipi{zB) ^ iPi{zb) 
S.i{zh)'^1{zh) = 



J+i 



(2^ + 1)!!' 

(2/-1)!! 
-^ 1 ' 



(18) 
(19) 



the dispersion relation (14) is fulfilled for any 
radius R of the sphere, and leads to the relation: 



/ 



^{^) 



= 1 



^ + 1 Sout 

giving discrete plasmon frequencies: 



^0,1 



I 



21 + 1 



(20) 



(21) 



which are real, in contrary to the exact solutions 
presented in figures 1 and 2 which arc obligatory 
complex. 

LUi{R) dependence resulting from the exact so- 
lution do not smoothly tend to the value lui{R -^ 
0) with decreasing R, as usually expected (e.g. [3], 
[5], [10]), but it grows up to Wi(i?min,i) value, as 
illustrated in figures 1. For R < Rmin,i there are 
no eigenvalues ili{R). This behavior of the il.i{R) 
dependence is mainly due to fast divergence of 
the ^i{zh) function entering the dispersion relation 
(14) for the arguments smaller than the range of 
variability of zh = zh{^i{R), R) parameters for 
successive I. 

When one includes the relaxation rate of the 
electron gas into the Drude model of the dielectric 



function, the plasmon frequency lui for given radius 
R of the sphere is relatively slightly red shifted, 
while to'/ experiences strong modification as illus- 
trated in figure 3 and 4 respectively for 7 = 1 eV 
([15]). 



the " low radius limit" approximation: ti>((i?niin^;) « 
cjp\/l/ {21 + 1). Even though the problem of the 
optical properties of metal sphere is at least as old 
as Mie theory [6], it seems, that the limitation for 
the smallest cluster still enabling the plasmon os- 
cillations has not been discussed previously. 



4. Conclusions 

By carefully studying the radius dependence of 
eigenmode problem of a sphere one can formulate 
several conclusions allowing for better understand- 
ing of surface plasmon features. In this paper we 
concentrate on studying the differences of surface 
plasmon features in the classical picture resulting 
from treating the radius dependence exactly, and 
the expectations from the widely applied approxi- 
mation of the so called " low radius limit" . We use 
the example of sodium sphere of plasma frequency 
Up = 5.6 eV , however the conclusions are quali- 
tatively valid for other simple free-electron metals. 
According to the non-approximated treatment the 
surface plasmons are always radiatively damped, 
even in the absence of coUisional process: eigenfre- 
quencies must be complex. The " low radius limit" 
leads to the real eigenfrequencies loq.i, which are ra- 
dius independent. From the exact calculations one 
can conclude, that the radius dependence of multi- 
polar plasmon frequencies is more subtle, than ex- 
pected. Our calculations show, that at larger polar- 
ity the uji{R) dependence docs not smoothly tend 
to the value wq,; = LOp\Jlj {21 + 1) of the vanishing 
size limit, as one could expect (e.g. [12], [13]). If the 
sphere is of radius R smaller than the characteristic 
radius Rminj ~ [I {21 + 1)] , there is no related 
eigenvalue ili {R) real nor complex. So the complex 
eigenfrequencies r2;(i?) = lui{R) + iu!i{R) can be 
attributed to the sphere starting from the radius 
R = -Rmin,; 7^ 0. The radii -Rmin,; for higher polari- 
ties / are not much smaller then the wavelength of 
the optical range (the anomalous dispersion range 
of alkalies) so the " low limit approximation" loses 
its validity. Our "numerical experiment" proves, 
that for the smallest particle radius -Rmin,; still pos- 
sessing an eigenfrequency in given polarity /, the 
plasmon oscillation frequencies can be well approx- 
imated by the corresponding value resulting from 
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